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Abstract 



QCD lattice simulations with 2 + 1 flavours typically start at rather 
large up-down and strange quark masses and extrapolate first the strange 
quark mass to its physical value and then the up-down quark mass. An 
alternative method of tuning the quark masses is discussed here in which 
the singlet quark mass is kept fixed, which ensures that the kaon always 
has mass less than the physical kaon mass. It can also take into account 
the different renormalisations (for singlet and non-singlet quark masses) 
occurring for non-chirally invariant lattice fermions and so allows a smooth 
extrapolation to the physical quark masses. This procedure enables a wide 
range of quark masses to be probed, including the case with a heavy up- 
down quark mass and light strange quark mass. Results show the correct 
order for the baryon octet and decuplet spectrum and an extrapolation to 
the physical pion mass gives mass values to within a few percent of their 
experimental values. 
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1 Introduction 



There has been a steady progression of lattice results from a quenched sea to a 
two-flavour and more recently 2 + 1 flavour sea in an attempt to provide a more 
complete and quantitative description of hadronic phenomena. (By 2 + 1 flavours 
we mean here two mass degenerate up-down, mf , quarks and one strange, , 
quark.) In this letter we discuss some ways of approaching in the mf-m^ plane 
the physical point {mf*,m'^*), where the natural starting point for these paths 
is an SU{?)) flavour symmetric point mf = = mfjJJ]. (The superscripts *, 
denote the physical point and flavour symmetric point respectively and ^ means 
the renormalised quantity.) The usual procedure is to estimate the physical 
strange quark mass and then try to keep it fixed, i.e. mf = constant, as the light 
quark mass is reduced to its physical value. However the problem is that the 
kaon mass is always larger than its physical value. We propose here instead to 
choose the path such that the singlet quark mass is kept fixed, 

m"* = 5(2mf + mf) = constant . (1) 

This procedure has the advantage that we can vary both quark masses over a wide 
range, and is thus particularly useful for strange quark physics. SU{?))f chiral 
perturbation theory should work well, because both the kaon and rj are lighter 
than their physical values along the entire trajectory. (They both approach their 
final mass values from below.) Since SU{2)) p chiral perturbation theory is thought 
to be valid for itik < 600 MeV, [1] , we should always be able to make use of chiral 
perturbation theory. If we extend our measurements beyond the symmetric point 
we can also investigate a world with heavy up-down quarks and a lighter strange 
quark. 

As a 'proof of concept' results given here show firstly the correct order for the 
baryon octet and decuplet mass spectrum and secondly an extrapolation to the 
physical pion mass yields results to within a few percent for the baryon masses. 



2 Extrapolating flavour singlet quantities 



Flavour singlet quantities are flat at a point on the SU (3) flavour symmetric line 
and hence allow simpler extrapolations to the physical point. This may be shown 
by considering small changes about a point on the flavour symmetric line. Let 
X5(mf,mf,mf) be a flavour singlet object i.e. Xs is invariant under the quark 
permutation symmetry between m, d and s. So Taylor expanding Xs about a 
point on the symmetric line where flavour 5't/(3) holds gives 



Xs(m«(°) + 5mf,m^(°) + 5mf,m«(°) + 5m 
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But on the symmetric line we have 
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sym 
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sym 
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sym 



(3) 



and on our chosen trajectory = constant, 

25mf + 5mf = , (4) 

which together imply that 

Xs(m«(0) + 5mf , m^(°) + 5mf, m^(°) + (5mf ) = xfly^ + 0{{5m^f) . (5) 

In other words, the effect at first order of changing the strange quark mass is 
cancelled by the change in the light quark mass, so we know that Xs must 
have a stationary point on the SU{3)f symmetric line. If we were making a 
quadratic extrapolation to the physical point, this eliminates a free parameter 
- we only need two parameters to make the quadratic extrapolation. This is 
particularly useful for quantities like the force scale, tq, where we do not have 
any theoretical input from chiral perturbation theory. Also the fact that Xs is 
flat at the symmetric point means that the extrapolated value cannot lie very 
far from the measured value. Other paths are not so fortunate. What we are 
doing is to keep the flavour singlet part of the quark matrix constant, while 
increasing the octet part (the changes in quark masses are proportional to Ag). 
Gluons are flavour-blind, they cannot couple directly to a flavour octet operator, 
they can only see flavour singlets - which can only occur in the square of the 
octet part of the mass. So everything about our gluon conflguration will vary 
quadratically with the distance from the symmetric point. So we are generating 
all our conflgurations closer to the physical point by taking the = constant 
line. 

Other potential advantages include: as mf \ mf* then mf iTif*, i.e. , the 
m^-mf splitting or increases to its physical value. The singlet quark mass is 
correct from the very beginning; numerically the simulation cost change should 
be moderate and the update algorithm is expected to equilibrate quickly along 
this path. 

For Xs we have several possibilities: for example the centre of mass squared 
of the meson octet, ^{m^ + 2m|f) or the centre of mass of the baryon octet or 
decuplet, Xn = ^(mAr+mE+ms) = 1.150 GeV, Xa = l{2mA+mn) = 1.379 GeV 
respectively or a gluonic quantity such as Xr = I/tq. 

We can check the above result, eq. ([5]), by considering leading order (LO) 
together with next to leading order (NLO) SU{3)f chiral perturbation theory, 
xPT. Note that LO xPT corresponds to linear terms in the expansion of Xs and 
so from eq. will be absent. Rather than using the full symmetric 1 + 1+1 results 
it is sufficient to just consider the 2 + 1 results. Let us first define Xg = 2B^m'^ 
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and furthermore set = ixi + 2Xs)/3 and Xk = {xi + Xs)/2. Then xi + Xv = 2x 
and Xi + '^Xk = 3% are constants on our trajectory and so 

Sxi + 6xrj = = 6xi + 26xk ■ (6) 
This means that any functions of the form 

fiXn) + fiXi) or 2h{xK) + Hxi) , (7) 

will have zero derivative at the symmetric point, so they are permitted as higher 
order corrections. Using the LO and NLO results from e.g. [2l|3lll] (where further 
details of the functions can be found) we have 

l{2ml + ml) = x + {fAx,) + Mxi)} 
i(mAr + ms+m=) = ttion + 2{aN + Pn + 3aN)x 

+{fNiXv) + fNiXi) + 2hN{XK) + hN^xi)} 
\{2mA+m^) = + 2(7a - 3o-a)x 

+{/a(x,) + /a to) + 2hAixK) + hAixi)} , (8) 

with 

U ix) = + PnX^ + 7- ln(x/Ax) , (9) 

and 

fsix) = SsX^+^sF{xh 

hsix) = CsX^+VsFixh: (10) 

(with S = N, A), a, . . ., 7] are combinations of the low energy constants. The 
NLO results are shown in curly brackets. These results are thus all in agreement 
with our previous discussion. Note also that on the = constant trajectory 
some of the low energy constants combine together leading to fits with fewer free 
parameters. 

We now have to relate the known physical point to the initial symmetric point. 
As discussed above, we expect Xg to be constant (in mf ) up to small corrections, 
so we will find it sufficient to consider only LO xPT. Then keeping constant 
means keeping {2mj^ + m^)/3 constant. For this path choice, we can now relate 
the known physical point to the initial symmetric point, 

l{2ml + ml) 

where S = N, A and r respectively. So simulations along the flavour symmetric 
line and using eq. ( ITTi) are sufficient to determine the initial point. This procedure 
is illustrated in Fig. [1] in the next section. 



mz. 

7r 
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(11) 



5 



If we now generalise to consider higher order terms from the above discussion 
lines of constant Xs are curved though they do still have to have the slope of —2 
at the point where they cross the mf = = mf^^ line. Now we have to specify 
more closely what we mean when we keep 2m|' + constant, as different scale 
choices give different paths. As shown above in eq. on our trajectory this is a 
higher order effect, on other choices of trajectory the difference between different 
scale definitions would be much more important. If we make different choices of 
the quantity we keep constant at the experimentally measured physical value, for 
example the choices discussed above we get slightly different trajectories. The 
different trajectories begin at slightly different points along the symmetric line. 
Initially they are all parallel with slope —2, but away from the symmetry line 
they can curve, and will all meet at the physical point. 



3 Clover fermions 

The above properties are general; we now apply them to 'clover' or 0(a)-improved 
Wilson fermions. There are now some additional points to consider when relating 
the bare quark mass (which is the input parameter) to the renormalised quark 
mass. The problem is that for fermions with no chiral symmetry the singlet, S, 
and non-singlet, NS, quark mass can renormalise differently which means that the 
relation to the bare quark masses and hence k, which is the adjustable simulation 
parameter, is more complicated |5j 

= Z^%mg + azm), (12) 
{q = I, s) with az = {Z^ — Z^)/Z^ and bare quark mass defined by 

"""^'^K^'^^) ' ^^^^ 

where Ksym;c is defined by the vanishing of the quark mass along the symmetric 
line, i.e. for 3 mass degenerate flavours. Now from LO xPT we have 

\{2{amKf + {am^f) oc |(1 + az)am , (14) 

and so the path am = constant remains as ^{2{amKy + {am^^y) constant. This 
translates to 



.(0) 



1^1 

ym 



where k^^J^ is the appropriate k on the SU{3) flavour symmetric line. As discussed 
previously higher order corrections have zero derivative at the flavour SU{3) 
symmetric point and so should be small. 
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0(a)-improvement also leads to a change in the coupling constant, gQ^gQ = 
(7q(1 + bgOm), [6l [7]. However for our trajectory as m is held constant then 
also remains constant. Also, [7], = [1 + {3dm + drnu{^))am]fn where dm, 
dm are improvement coefficients and m(^) = 3 — 4^ + 2^^ with ^ = mi/mf^^. As 
m(^) has a minimum at ^ = 1, it is again flat at the symmetry point, as expected. 
We have estimated that with the parameters used below, the correction term in 
j^R j^gjxiains very small and so we shall ignore this here. 

The particular clover action used here consists of the tree level Symanzik 
improved gluon action together with a mild 'stout' smeared fermion action. Fur- 
ther details, [H |9], and a determination of the non-perturbative, NP, coefficient 
used for the clover term are described in [9]. Simulations have been performed 
using the Hybrid Monte Carlo, HMC, algorithm with mass preconditioning for 
2 mass-degenerate flavours and the rational HMC, [TU] for the 1-fiavour. Two 
programmes were used, a Fortran programme, [TT], and also the Chroma pro- 
gramme, [12]. All the runs described here are on 24^ x 48 lattices at /3 = 5.50, 
Csw = 2.65, [9] (we describe the determination of the scale later). We ffist need 
to determine the symmetric point, using eq. ffTTl) . 

A series of runs on the SU (3) flavour symmetric line gives an estimate of our 
starting value n^^Jm- Some are shown in Fig. [1] (where we plot (2m|- — m^)/X| 
against m^/X| using the scales Xs with S = N, A and r) as the points lying on 
the y = X = m^/X|, black dashed line. We seek the point on the symmetric line 
where eq. fill I) holds. This gives an estimate for We find k^^J^ = 0.12090 

which we shall take as our starting value. 

From this we can now, given a ni, find the corresponding Kg using eq. f|T5|) . 
After some experimentation we chose the ki, Kg values given in Table [TJ Note 





Kg 




0.12083 


0.12104 


mi > nig 


0.12090 


0.12090 


mi = mg 


0.12095 


0.12080 


mi < mg 


0.12100 


0.12070 


mi < mg 


0.12104 


0.12062 


mi < mg 



Table 1: {Ki,Kg) values simulated on 24^ x 48 lattices. 

that it is possible to choose k/, Kg values (here (0.12083,0.12104)) such that 
mi > mg. In this strange world, as previously mentioned, we would expect to 
see an inversion of the particle spectrum, with for example the nucleon being the 
heaviest octet particle. 

These result^ are also shown in Fig. [TJ Each data set comprises ~ 0(2000) 
trajectories. (Note that our lowest pion mass has m^L ~ 3.4 and may be showing 

^Preliminary results were given in |13j . 
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m '/X ' 

Figure 1: (2m|^ — m^)/X| (y-axis) against m'^/X'g (x-axis) using S = N (red circles) 
5 = A (green squares) and S = r (blue diamonds). The SU{3) flavour symmetric line 
{y = x) is the dashed line. (For convenience the results for S = r have been divided by 
a factor of 20.) The experimental points using the three singlet quantities, Xs, S = N, 
A, r are shown as stars. The solid lines are fits using eq. (|16p . 



some sign of finite size effects.) Also shown is a fit to constant 2mj^ + by 

^^-Cs-2<. (16) 

The numerically simulated points all lie (approximately) on the line of constant 
2m|; + m^. There is also consistency between the various singlet quantities used. 

To determine the scale we again use the constancy of Xs- In Fig. [2] we show 
aXs against (am,r)^ for S = N, A, r together with constant fits. (Although 
the physical limit is also shown, as we are making a constant fit, this is not 
important here.) Then from aXs = constant we can determine the scale giving 
a = 0.083 fm, 0.084 fm for S = N, A respectively. We use in future the S = N or 
a = 0.083 fm value. This means that the box size, L ~ 2fm. While for X = r the 
numerical results are the flattest, the physical value is less well known. However 
reversing the argument and using a = 0.083 fm gives tq = 0.50 fm. 

As is apparent from Fig. [H we have slightly underestimated k^^J^ (reflected 
by the fact that ^cs from eq. ( |T6i) is not quite equal to |(2m|' + m^)/X||*, cf 
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Figure 2: aXs against {am.,^)'^ for S = N (red circles) = A (green squares) and 
S = r (blue diamonds). Constant fits are also shown (dashed lines). 



eq. ( fTTl) ). We can estimate the significance of tliis on the kaon mass, by taking the 
value of {2m'j^—m'^)/Xff at the physical pion mass, (m^/X^)|* in the figure. This 
gives ttlk ~ 509 MeV, a ~ 3% discrepancy when compared to the experimental 
kaon mass of ~ 494 MeV. 



4 Hadron spectrum 

The octet (and decuplet) baryon masses are degenerate at the SU{3) flavour 
symmetric point and then fan out. As a first example we now give some mass 
results in Fig. |3]for the baryon octet: A^, A, E, S, together with a constrained fit 
from LO xPT, e.g. [3] about the SU{3) flavour symmetric point, 

rriN = AN + 2{aN + (3N){ml-ml\^^J.J 

rriY. = An + (oat - 2/3Ar)(m^ - ^llfjm) 

tua = Aat - (ftTv - 2/3Ar)(m^ - m^l^L) 

ms = An - 3aN{ml - ml\fjj , (17) 
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0.05 0.1 0.15 0.2 



Figure 3: The octet baryon masses O where O = A^, A, S and H using the scale 
S = N, upper left to lower right respectively, together with the fit from eq. (jl7p . The 
common SU{3) p symmetric quark mass value is shown in red. The experimental values 
are shown with stars. 



where A]^ = itlqn + '^{q-n + Pn + 3aN)Tn, aN, Pn are the three fit parameters 
(in the figure we have also normaUsed the masses with Xn). Similarly in Fig. H] 
results are given for the baryon decuplet: A, E*, S*, f2, again together with the 
constrained fit from LO xPT, e.g. [1], 

m-s* = Aa 

m^, = Aa - 2-fA{ml - mllfjj 

mn = Aa - A^Aiml - ml\%l) , (18) 

where A a = w^oA + 2 (7a — 3(rA)m, 7a are the two fit parameters. 

Both these figures illustrate the 'proof of concept' of the method described 
here. The correct ordering of the particle spectrum has been achieved (including 
the anti-ordering behind the symmetric point). The masses (using the scale 
determined by 5* = N) are given in Table [2l The results are already within a few 
percent of their experimental values, and clearly show that higher orders of xPT 
are small. We can expect an improvement for results closer to the physical pion 
mass. 
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Figure 4: The decuplet baryon masses, D, where D = A, S*, H* and 0, (also using 
the scale S = N), together with the fit from eq. (jlSp . The common SU{3)f symmetric 
quark mass value is shown in red. The experimental values are shown with stars. 



Finally to illustrate better the 'fanning' out of the results we show in Fig. [5] 
mass splittings for the baryon octet: m]\f — m^, mA — rriY:, — mj^ together with 
the constrained fit from eq. (fTTl) . Considering mass splittings has the advantage 
that the results can be obtained directly from the ratio of correlation functions, 
which leads to a significant reduction in the error. 



5 Conclusions 

We have suggested here that the simplest way of approaching the physical point 
in lattice simulations is to hold the singlet quark mass fixed from a point on 
the SU{3) flavour symmetric line. This has been shown both theoretically and 
numerically (using an NP 0(a)-improved 2 + 1 flavour clover action) to lead 
to very smooth results in the extrapolation of singlet quantities to the physical 
pion mass. Exploratory results for the hadron mass spectrum give masses in 
the correct order (including inversion when mi > rris i.e. we can simulate a 
strange world where, for example, the nucleon can decay). Furthermore the 
extrapolated masses for both the baryon octet and decuplet are within a few 
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ID, AT 


939 


998(16) 
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rriA 


1116 


1133(14) 
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ms 
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ms 


1318 


1318(15) 
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mA 


1232 


1303(16) 


5% 


rriY:* 


1385 


1402(05) 


1% 


ms* 


1533 


1501(12) 


2% 


mn 


1673 


1599(25) 


4% 



Table 2: Masses for the baryon octet and decuplet. 



percent of their experimental values. To improve the situation further clearly we 
need simulations closer to the physical pion mass. At present in our simulations 
the pion mass decreases from ~ 450 MeV to ~ 335 MeV while the kaon mass 
increases from ~ 400 MeV to ~ 450 MeV, (experimentally m^r = 138 MeV, rriK = 
494 MeV). Also to improve the accuracy of the approach to the physical pion 
mass another line of constant singlet mass should be found, which will allow 
interpolation/extrapolation around the physical pion mass. Finally we note that 
even LO xPT seems to be working very well around the flavour symmetric point. 
Further results will be published elsewhere, [Ti] . 
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